Radiation-Driven Outflows in Active Galactic 

Nuclei 



Daniel Proga* and Ryuichi Kurosawa*'* 

* Department of Physics and Astronomy, University of Nevada Las Vegas, Box 454002, 

4505 Maryland Pkwy, Las Vegas, NV 891541-4002 
' Current Address: Department of Astronomy, Cornell University, Ithaca, NY 14853-6801 

Abstract. We review the results from multi-dimensional, time-dependent simulations of gas dy- 
namics in AGN. We will focus on two types of outflows powered by radiation emitted from the 
AGN central engine: (i) outflows driven from the innermost part of an accretion disk and (2) out- 
flows driven from a large-scale inflow that is likely the main supplier of material to the central 
engine. We discuss the relevance of both types of outflows to the so-called AGN feedback problem. 
However, the AGN feedback should not be considered separately from the AGN physics. Therefore, 
we also discuss the issue whether the properties of the same outflows are consistent with the gas 
properties in broad- and narrow-line regions. 
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INTRODUCTION 

The dynamics of the gas and dust in narrow line regions (NLRs) and broad line regions 
(BLRs) in active galactic nuclei (AGNs) is driven by gravity but radiation driving 
can also be important even for sub-Eddington sources (e.g., yj]). Radiation driving 
can be due to radiation pressure, radiation heating, or both. The radiation force can 
overcome gravity for sub-Eddington sources when the gas/dust opacity is higher than 
the electron scattering. The latter is usually used to define the Eddington luminosity, 
^Edd- The gas opacity can be enhanced by the scattering of photons by UV spectral 
lines. Radiation pressure on spectral lines (line force) can be significant, provided that 
the gas is moderately ionized and can interact with the UV continuum through very 
many UV line transitions. For highly ionized gas, line force is negligible because of a 
lower concentration of ions capable of providing UV line opacity. In the case of highly 
or fully ionized gas, an outflow can still be produced if the gas heating is efficient enough 
for the thermal energy to exceed the gravitational energy. 

AGNs with their broad spectral energy distributions (SEDs), are systems where both 
line driving and radiation heating, in particular X-ray heating, can operate. In fact, a wind 
driven from an accretion disk by line force is the most promising hydrodynamics (HD) 
scenario for outflows in AGN, especially high-luminosity quasars. In this scenario, a 
wind is launched from the disk by the local disk radiation at radii where the disk radiation 
is mostly emitted in the UV ([|20; Q3U ) . Such a wind is continuous and has a mass-loss 
rate and velocity that are capable of explaining the blueshifted absorption lines observed 
in many AGNs, if the ionization state is suitable (e.g., [H]; 0]; This wind scenario 



has a desirable feature, i.e., for the wind motive power it relies on radiation, which is an 
observable quantity. We note that line-driving can account not only for gross properties 
of some AGN outflows but also for specific spectral features observed in some quasars 
id])- However, not all AGN outflows can be explained by line driving because of 
too low a luminosity, too high an ionization state, or both (e.g., U; lioll ). Therefore, 
other mechanisms such as thermal and magnetic driving are likely also important. 

Theoretical models predict that X-ray heating can have profound effects on the gas 
dynamics in disks. Since X-rays tend to heat low-density gas to a temperature Tq ~ 
10 7 K, with which matter in an accretion disk is expected to either puff up and form 
a static corona or produce a thermal wind, depending on whether the thermal velocity 
exceeds the local escape velocity, v esc (e-g-, [El]; CLll; 10]; tfHl])- 

These studies demonstrate that radiation liberated by an accreting disk can drive 
a powerful outflow from this disk. Therefore one can ask: Can AGN radiation drive 
an outflow from anywhere other than the disk? If yes, what are the properties of this 
outflow? In particular, can it explain the properties of BLR and NLR? Can this outflow 
be a part of AGN feedback? To answer these questions, we recently applied and extended 
our techniques developed to study disk outflows to model fluid dynamics on scales 
comparable to the BH's gravitational radius and also on larger scales ( rfl5ll : lfl6ll : lfl7ll : 
[fl8|l : [fl9|l : ||20(l ). In the following, we briefly summarize the main results of our recent 
work. 



SIMULATIONS OF OUTFLOWS FROM INFLOWS 



In H15H . we calculated a series of models for non-rotating flows that are under the in- 
fluence of super massive BH gravity and radiation from an accretion disk surrounding 
the BH. Generally, we used the numerical methods developed by B4Q. Our numerical ap- 
proach allows for the self-consistent determination of whether the flow is gravitationally 
captured by the BH or driven away by thermal expansion or radiation pressure. 

To compute the structure and evolution of inflows/outflows, we solve the equations of 
HD: 

Dp 

^ +P V.v = 0, (!) 
p — = -V/>-pVcD + pF rad , (2) 

PTz( e ~l) =- J PV-v + pif, (3) 



Dt VP 

where p is the mass density, P is the gas pressure, v is the velocity, e is the internal 
energy density, Jzf is the net cooling rate, 4> is the gravitational potential, and F rad is the 
total radiation force per unit mass. 

In our previous simulations, the total radiation force has two components: one is 
due to electron scattering, F rad ' e and the other is due to lines F 18 * 1 ' 1 . The latter can be 



approximated a modified version of the method developed by |21|1. The line force at a 
point defined by the position vector r is 



F"*'(r) = lujt) ( h °^ da ) (4) 



where / is the frequency-integrated continuum intensity in the direction defined by the 
unit vector n, and £1 is the solid angle subtended by the disk and corona at the point. The 
term in brackets is the electron-scattering radiation force, o e is the mass-scattering co- 
efficient for free electrons, and M(t) is the force multiplier - the numerical factor which 
parametrizes how much spectral lines increase the scattering coefficient. In the Sobolev 
approximation, M(t) is a function of the optical depth parameter t = a e pvth/\dvi/dl\, 
where v t h is the thermal velocity, and dvi/dl is the velocity gradient along the line of 
sight, h. The force multiplier approaches zero when all lines become optically thick 
(i.e., lim r ^ 00 M(?) = 0) and some finite maximum value (M max ) when all lines become 
optically thin (i.e., lim/_>o = ^max)- The maximum value of the force multiplier is a 
function of physical parameters of the wind and radiation field, and can be parametrized 
by the photoionization parameter, £ . Several studies have shown that M max is roughly 
a few thousand for gas ionized by a weak or moderate radiation field (e.g., Q2JJ] ; |22|1 : 
lHH]; [0]). As the radiation field becomes stronger and the gas becomes more ionized 
the force multiplier decreases asymptotically to zero (i.e., lim^ M max = a few 10 3 and 
lim^^Mmax = 0). M max allows us to set a lower limit of the luminosity at which a sys- 
tem can have a "line-driven" wind, [i.e., L > Lg/{ \ +M max )]- We note that using the 
Castor et al.'s method 112 ill , one can model large scale properties of relatively smooth 
line-driven winds very well, but cannot capture small scale effects, in particular the so- 
called line driven instability (e.g., [25]; (20]; §2% }. 

We adopt an adiabatic equation of state P = (y— l)e, and consider models with 
the adiabatic index, 7= 5/3. We performed simulations using the Newtonian potential 
4> due to the central BH (the general relativity effects can be neglected because we 
consider flow dynamics relatively far from the BH). The simulations were performed 
in spherical polar coordinates assuming axial symmetry about the rotational axis of the 
accretion disk. 



RESULTS 

For a 10 8 M BH with an accretion luminosity of 0.6 of L^dd, we found that a non- 
rotating flow settles quickly into a steady state and has two components (1) an equatorial 
inflow and (2) a bipolar inflow/outflow with the outflow leaving the system along the 
pole. The first component is a realization of Bondi-like accretion flow. The second 
component is an example of a non-radial accretion flow becoming an outflow once it 
is pushed close to the rotational axis of the disk where thermal expansion and radiation 
pressure can accelerate the flow outward. The main result of these simplified calculations 
is that the existence of the above two flow components is robust yet their properties 
are sensitive to the geometry, SED of the radiation field, and the outer boundaries. In 
particular, the outflow power and the degree of collimation are higher for the model with 
radiation dominated by UV/disk emission than for the model with radiation dominated 
by X-ray/central engine emission. This sensitivity is related to the fact that thermal 
expansion drives a weaker and wider outflow, compared to the radiation pressure. 




FIGURE 1. Three-dimensional hydrodynamical simulations of outflow formation via redirection of 
accreting gas by the strong radiation from an accretion disk around a super massive black hole with its 
mass M BH = 10 s . VI.. The inf ailing gas is weakly rotating (sub-Keplerian), and the Eddington ratio of 
the system is 0.6. The volume rendering representation of the density distributions is shown. The outflow 
morphology is bi-conical, but the flow contains relatively cold and dense cloud-like structures which 
resembles those observed in the NRLs of Seyfert galaxies. The figure is from ifTill . 



In iflfjll . we explored effects of gas rotation. As expected, rotation changes the geom- 
etry of the flow because the centrifugal force prevents gas from reaching the rotational 
axis (see fig. 1 in 03]). 

This, in turn, reduces the mass outflow rate because less gas is 
pushed toward the polar region. We also found that rotation can lead to fragmentation 
and time variability of the outflow. As the flow fragments, cold and dense clouds form 
(figs. 4 and 12 in £3]). 

We have also started studying 3-D effects on gas dynamics (e.g., lfl7ll : lll8ll ). In 1 171. 



we considered effects of radiation due to a precessing accretion disk on a spherical cloud 



of gas around the disk On the other hand, in [18]), we recalculated some models from 
papers ||15|l and dl6fl in full 3-D. Our 3-D simulations of a nonrotating gas show small 
yet noticeable nonaxisymmetric small-scale features inside the outflow. However, the 
outflow as a whole and the inflow do not seem to suffer from any large-scale instability. 
In the rotating case, the nonaxisymmetric features are very prominent, especially in 
the outflow which consists of many cold dense clouds entrained in a smoother hot flow 
(e.g., see figs. [TJ and [2]) . The 3-D outflow is nonaxisymmetric due to the shear and thermal 
instabilities. Effects of gas rotations are similar in 2-D and 3-D. In particular, gas rotation 
increases the outflow thermal energy flux, but reduces the outflow mass and kinetic 



Sample movies of the 3-D simulations can be found at |http : / /www . physics . unlv . edu/~rk/ research/agn_3cL3~'' 
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FIGURE 2. A scatter plot of the radial velocity (v,-) verses gas temperature (T) for the model shown in 
fig-El A negative value of v r indicates an inflow. To avoid overcrowding, only the points on <p = plane 
are shown. A large fraction of gas is in outflow motion (v, > 0) for the models with rotation and a wide 
range of the temperature is associated with the outflowing gas. The figure is from lfl8ll . 



energy fluxes. In addition, rotation leads to time variability and fragmentation of the 
outflow in the radial and latitudinal directions. The collimation of the outflow is reduced 
in the models with gas rotation. The main different effect of rotation in 3-D compared to 
2-D is that the time variability in the mass and energy fluxes is reduced in the 3-D case 
because of the outflow fragmentation in the azimuthal direction. 

Our recent simulations, demonstrate that AGNs can have a substantial outflow orig- 
inating from the infalling gas. Such an outflow can reduce the rate at which matter is 
supplied to the central region of AGN because its mass loss rate can be significantly 
higher than the mass inflow rate at small radii (see fig. [3]). For example, as little as 10% 
of the inflow at large radii can reach small radii because 90% of the inflow is turned 
into an outflow. In general, the kinetic power dominates the thermal power at all radii in 
both models with and without gas rotation; however, the thermal power contribution is 
non-negligible in the models with gas rotation (fig. [3]). The X-ray heating becomes more 
effective in the rotating gas environment because the flows become more non-radial and 
will be subjected to the direct exposure to the strong radiation from central continuum 
radiation source. 

In [fl9|l . the radiation-driven AGN outflow model of [T^\ is extended by relaxing the 
assumption of a constant accretion luminosity. This allows us to determine the accretion 
luminosity consistently with the mass accretion rate at the inner boundary, and conse- 
quently the two quantities are coupled through the radiation field. This is an improve- 
ment toward a more comprehensive self-consistent hydrodynamical models with radia- 
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FIGURE 3. Mass and energy fluxes as a function of radius for the model shown in fig. Q] The panel 
is subdivided into two parts: top (mass flux) and bottom (energy flux). In the mass flux plot, the inflow 
(dashed line; M; n ), outflow (solid line; M Q ) and net (dotted lines; M net ) mass fluxes, are separately plotted, 
as a function of radius. The absolute values of M; n and M net are plotted here since they are negative at 
all radii. The length scale is in units of the inner disk radius (r' — r/r* where r* = 8.8 x 10 13 cm). In the 
energy flux plots, the kinetic energy (solid line) and the thermal energy (dotted line) fluxes are shown. The 
figure is from IU8I1 . 



tive feedback. For the models with high temperature gas at large radii (~ 10 pc) and high 
luminosities, we find a strong correlation between the mass-outflow rate (M ut) and the 
luminosity (L). The power law index (q) describing the M out -L relation is q = 2.0 (±0.1), 
which is very similar to that for radiation-driven stellar and disk wind models. More sur- 
prisingly, for high density at large radii, we find steady state solutions with the accretion 
luminosity exceeding the Eddington limit. The super-Eddington accretion proceeds in 
the equatorial region and is possible because the radiation flux from the disk is signifi- 
cantly reduced in the equatorial direction due to the geometrical foreshortening effecll 
(see fig-IU). 

Based on a set of axisymmetric simulations presented in Ufa] , we analyzed the energy, 
momentum, and mass feedback efficiencies due to radiation from AGN (see [20]). We 
find that even for the strongest outflow, the ratio between the outflow kinetic power and 
the radiation power is very low (~ 10~ 4 at the peak; see fig. [5]). One of the reasons for 
this relatively low efficiency of the large scale outflows is that for large densities at large 
radii, accretion proceeds in the equatorial region, and the radiation flux from the disc 
is significantly reduced the geometrical foreshortening effect, as also mentioned above. 



Sample movies of the 2-D simulations can be found at |http : / /www . physics . unlv . edu/~rk/ research/agn2d_i"u 
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FIGURE 4. Examples of the density and velocity maps from a set of two-dimensional simulations 
(without the temperature constrained at the outer boundary) presented in 11911 . The density (in logarithmic 
scale) is over-plotted with the directions of poloidal velocity as arrows. The figures are placed in order 
of increasing density at the outer boundaries (p G ) from the left to right. The self-consistently determined 
accretion luminosities (in the unit of Eddington luminosity) are F = 0.32,0.71 and 4.3, from the left 
to right panels. For the low density (p G ) and low accretion luminosity model (left panel), the outflow 
is very narrow (0° < 9 < 30°), and the inflow is very wide (30° < 9 < 90°). As the density (p G ) and 
accretion luminosity increase (middle panel), the outflow becomes wider (0° < 9 < 50°), whereas the 
inflow becomes narrower (50° < 9 < 90°). For the very high density and accretion model (right panel), 
the outflow occurs over a very wide range of the polar angle (0° < 9 < 85°), and the accretion region is 
now confined to a thin equatorial wedge (the disk-wind-like solution). The figures are from Q20I1 . 



The coupling between the radiation and matter becomes less efficient once the inflow- 
outflow morphology becomes the "disk wind like" (the right panel in fig. H]) because of 
the "mismatch" between the direction in which the disk radiation peaks (in the polar 
direction) and the direction of the matter inflow (in the equatorial direction). 

Compared to the energy (thermal only) feedback efficiencies ~ 0.05 required in the 
recent cosmological and galaxy mergers simulations (e.g., ([281] : |2§]; [f3Qh : |3l|l : 
our thermal energy feedback efficiency at the peak value is about 5 x 10 3 times smaller. 
Our total and kinetic energy efficiencies are about 5 x 10 2 times smaller than the value 
required in cosmological simulations. These large discrepancies might indicate a few 
things. For example, our models are missing important elements. In particular, we do 
not include effects of dust which could make the outflows much stronger. 

On the other hand, it is also possible that the AGN feedback may not be as effective as 
one might have had expected. Instead, other forms of feedback may be more significant 
than the AGN feedback via radiation on scales between 10~ 2 and a few parsecs. For 
example, the feedback via supernovae, star formation processes, the strong stellar wind 
from massive stars, and strong accretion disk winds (discussed in section 1) or jets from 
AGN may play more important roles. The last two forms of the feedback require a 
proper treatmentment of magnetic field, as they may carry a significant fraction of the 
total outward flux in energy and momentum. 

Finally, the AGN feedback efficiency may be indeed low, and the AGNs take a long 
time to influence their environment. We note that in our models the AGNs do not shut off 
the mass supply completely even at very high luminosities (T > 1). This indicates that 
the AGNs can operate on a very long time scale over which their affects can accumulate, 
and eventually become significant. 
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FIGURE 5. The efficiencies of converting the black hole accretion luminosity L a to the rate of energy 
deposition to the surrounding gas plotted as a function of the Eddington ratio (F). The values are based on 
the 2-D simulations similar to those shown in fig.|U The kinetic energy feedback efficiency (circles), 
the thermal energy feedback efficiency e t h (squares) and the total energy feedback efficiency e t = + £th 
(triangles) are shown separately. The maximum total energy feedback efficiency is ~ 10~ 4 . For the models 
with relatively low Eddington ratio (F < 0.4), the thermal feedback is more efficient than the kinetic 
feedback (e t h > £k)- For the models with relatively high Eddington ratio (r > 0.6), the kinetic feedback 
is more efficient than the thermal feedback by a factor of ~ 10 to ~ 100. The model with F = 0.2 does 
not form an outflow, and the vertical line (dashed) at F = 0.2 indicates an approximate F value below 
which no outflow forms. The flattening of the efficiencies beyond T « 1 is caused by the transition of the 
inflow-outflow morphology to a "disk wind like" configuration for the larger F models (c.f., fig. @). The 
figure is from l20ll . 

CONCLUSIONS 

The main conclusions from our recent simulations of radiation driven large scale out- 
flows in AGN include: 

• A significant fraction of the inflowing matter can be expelled by radiation heating 
and pressure. 

• The non-rotating flow settles into a steady inflow/outflow solution. Rotation of the 
flow and large optical depth lead to time variability in the inflow/outflow solution. 

• In the time variable flows, dense clouds form. This could be applicable to NLR in 
AGN. 

• The large outflows are multi-temperature/phase and could explain warm absorbers 
in AGN. 

• The mass supply rate does not appear to be limited by the AGN radiation. 

• The large-scale outflows are efficient in removing matter but inefficient in carrying 
out energy. 



• Radiation drives much more powerful outflows from small scales than from large 
scales. 
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